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Abstract Multiple zeta values (MZVs) are under intense investigation in three arenas - 
knot theory, number theory, and quantum field theory - which unite in Kreimer's proposal 
that field theory assigns MZVs to positive knots, via Feynman diagrams whose momentum 
fiow is encoded by link diagrams. Two challenging problems are posed by this nexus of 
knot /number/field theory: enumeration of positive knots, and enumeration of irreducible 
MZVs. Both were recently tackled by Broadhurst and Kreimer (BK). Here we report 
large-scale analytical and numerical computations that test, with considerable severity, 
the BK conjecture that the number, Dn,k, of irreducible MZVs of weight n and depth k, 

is generated by nn>3nfc>i(l - x'^y'')^"-'' = 1 - + -^rr^i^^^zj) , which is here shown 
to be consistent with all shuffle identities for the corresponding iterated integrals, up to 
weights n = 44, 37, 42, 27, at depths k = 2, 3, 4, 5, respectively, entailing computation at 
the petashuffle level. We recount the field-theoretic discoveries of MZVs, in counterterms, 
and of Euler sums, from more general Feynman diagrams, that led to this success. 



^) email: D.Broadhurst@open.ac.uk 



1. Introduction 



In a recent Physics Letter (hereafter BK), Dirk Kreimer and the author made a 
guess, informed by field theory, at the number, Dn ^, of irreducible multiple zeta values 
(MZVs) [@] of weight n and depth k that enter a minimal Q-basis, to which all other MZVs 
may be reduced, as rational combinations of elements of the basis set, and their products. 
MZVs, and their extension to alternating Euler sums, concerned BK, as practitioners of 
perturbative quantum field theory (pQFT), because of Kreimer's connection |, ||, of 
MZVs with positive knots 0, via the counterterms of pQFT. This connection is strongly 
supported by BK's joint |lO|, |l2l and separate ||13|, |14|, [T^ calculations of 



multi-loop Feynman diagrams, and by those of others workers [IS, 19, 20, 21 



At first parallel ||22|, ^ to, and recently intertwining p5| , p6| , p7t with, this flurry 



of field-theoretic activity, is an equally dramatic increase of understanding of MZVs and 
of alternating Euler sums, as richly structured number-theoretic objects, in their own 
right, irrespective of any field or knot theory. This is a rapidly progressing subject, where 



little is yet rigorously proven, but much is conjectured, on the basis of extensive P, K 



analytical and numerical computations. Several conjectures command wide @, ^, 

PB| , support. Yet the answer to a very obvious question - how many MZVs are left 

undetermined by the relations between MZVs? - appeared to be difficult, even to guess. 

As in 0, we approach the question in three stages. First, we consider Euler sums, of 
weight J2j Sj and depth /c, which are fc-fold nested sums of the form [jl], ^ 

k ^"'j 

C(si,...,Sfc;cri,...,crfc) = 

nj>nj+i>0 j=l 

with signs aj = ±1, positive integer exponents Sj, and aiSi 7^ 1, to prevent divergence 



of the outermost sum. As in p6[, we combine the strings of exponents and signs into a 
single depth-length argument string, with Sj in the jth position when aj = +1, and Sj in 
the jth position when aj = —1. The first question is then: 

Ql What is the number, En,k, of Euler sums of weight n and depth k, in a 
minimal Q-basis for reducing all Euler sums to basic Euler sums? 

Next, we define MZVs 0] as non- alternating Euler sums, with aj = 1, and ask: 

Q2 What is the number, Mn^k, of Euler sums of weight n and depth k, in a 
minimal Q-basis for reducing all MZVs to basic Euler sums? 

Finally comes the most natural, yet most difficult question: 

Q3 What is the number, Dn^k, of MZVs of weight n and depth k, in a 
minimal Q-basis for reducing all MZVs to basic MZVs? 

Section 2 gives the conjectured answers and focuses attention on the most difficult 
question, Q3, whose answer followed discoveries fl], ^ |ll| of how Feynman diagrams 
assign MZVs to knots. Section 3 characterizes the two types of identity that are believed 
to achieve the conjectured reductions: depth-length shuffies, and weight-length shuffies. 
Section 4 presents the results of (very) large-scale computational tests, both analytical 
and numerical. Section 5 presents a brief bestiary of MZVs and Euler sums, obtained from 
Feynman diagrams that played vital roles in the inference of the conjectured enumeration. 
A conclusion appears, with several open questions, in Section 6. 



1 



2. Conjectured enumerations 



It is conjectured that the answers to the three questions, above, are generated by 



n>3fc>l X X 1 

nn(i--v)"- = 1-1^' 

n>3fc>l 



The solution to (||) is that given by the author in P], in terms of what was there 
dubbed "Euler's triangle", with the symmetric entries: 

T(a, b) = ^ Y. M b/d) ' (4) 

+ " d\a,b 

where P{a,b) = (a + b)\/{a\b\) are the entries in Pascal's triangle, the sum is over all 
positive integers d that divide both a and b, and a Mobius transformation is effected by 

{1 when d = 1 

when d is divisible by the square of a prime (5) 

(—I)*'' when d is the product of k distinct primes 

which is the Mobius function. When n and k have the same parity, and n > k, one 
obtains 

^n,fc = T(^,A;). (6) 

With the exception of In 2 and tt^ (which act as seeds) all elements of the basis are thereby 
conjecturally enumerated. 

The solution to (0) is that given by BK. It too involves Euler's triangle. When n and 
k have the same parity, and n > 3k, one obtains 0] 

Mr,,k = T{^,k). (7) 

With the exception of vr^ and C(3) (which act as seeds) all elements of the basis are thereby 
conjecturally enumerated. 

Conjecture (|]) appears in the final version of BK, though no solution to it was 
given there. Here we motivate the Ansatz, and develop the generators for irreducibles of 
specific depths, up to A; = 6. Observe, first, that the order y terms in (|l]) merely assert 
the irreducibility of the 2-braid torus knot-numbers 0, C{2n + 1) for n > 0. The claim 
of [0, there supported by a year of effort and 10^ CPUhours of testing, is that all else is 
generated by including the factor of 1/(1 — xy) on the RHS, which gives (H), by Mobius 
inversion. It then took little time for BK to discover a RHS in (^ that was consistent 
with all known data. To celebrate the absence of In 2 from the reduction of both MZVs 
and counterterms, we omitted the factor of 1/(1 — xy), obtained (^, and performed CPU- 
intensive tests on the resulting prediction (|^), which emerged unscathed. It is important 
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to remember the question thus answered: how many Euler sums in the basis for MZVs? 
It is ironical that the most complex of the three questions admits of the simplest answer. 
Consider weight-12 depth-2 irreducibles, with Mi2,2 = 2. The two basis elements may 
be taken as C(9,3) and C(9,3). The latter is not an MZV, but should be included, since 
it was shown in that C(4, 4,2,2) reduces to a combination of C(9,3) with MZVs of 
depths k < 4. Moreover, at weight 15 and depth 5, we found that C(6,3,2,2,2) reduces 
to a combination of 14C(9,3,3) - 3C(7,5,3) with MZVs of depths k < 5. We call this 
phenomenon "pushdown" . The simplicity of relies on its existence and regularity. 

The final step, to the most complex of the three Ansatze, in (^, was taken more 
falteringly^ . The question - how many MZVs in the basis for MZVs? - is the most 
natural to ask, in knot /number/field theory One now suspects that its answer will 
be the most complex of the three, since one is denied the use of non-MZV terms, like 
C(9,3), in the reduction of MZVs. If pushdown makes (H) so simple, then its exclusion 
is expected to complicate (0). The final recipe of BK was the simplest possible, in these 
curious circumstances: to add a term that encapsulates what was already known about 
pushdown from depth 4 to depth 2, starting at weight 12. The additional term was 
constructed, with frank naivety, as follows. First, BK knew that the terms of order 
were given by (i_^4)(i_^6) ; merely by supposing the correctness of the conjectured [0, ^ 



enumeration of irreducible depth-2 MZVs. The factor of (1 — y'^) was then inserted on 
the grounds that pushdown to depth k = 2 originated from depth k = 4. This pushdown 
demonstrably 0, [Tl| occurs at weight 12, where the 'rule of three', -D2n,2 = L^^J) ^^r 
irreducible depth-2 MZVs, in [^, first differs from the 'rule of two', M2„.2 = L^T^J' f*^^ 



3-braid knots, in |jTT|. It is supposed that pushdown to depth k = 2 occurs for greater 



even weights, at a rate that is metered [|lT| by comparing [y] with |2J . 

It should be stressed that recipe (|^) is the simplest^ one that BK could concoct, on 
the basis of depth-2 results. Its genesis does not guarantee any success at depths k > 2. 
It is hoped that the reader will share some the author's amazement at the formidable 
success, below, of the highly specific predictions at depths k = 3,4, 5. To extract these 
predictions: take the logarithm of each side of (|^); set a = x"^ and b = x^y; expand in 
powers of 6; determine the generator for D^k^2n,k by Mobius inversion. The result is 

D,{a) = J2 ^3.+2n,fca" = Y: ^L./.ia") , (8) 

n>0 d\k 



where 



t^i (l-a2)(l-a3) 



TV 



(9) 



It follows that Dk{a) is a ratio of polynomials, and that its singularities occur exclusively 
at roots of unity. In particular: 

Di{a) = j^, (10) 
1 — a 

D2{a) = ^ -, (11) 

(1 — a)[l — a-^) 



preliminary version of BK got it wrong, by failing to model on the product forms 
^The aforementioned false start was more complicated, as it was not based on a product generator. 
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= J. \r, 3r 12 

(1 — a)\ \ — a j(l — a-^) 

. , _ 1 + 2a^ + + + 2a^ + - 

^'^""^ " (l-a)(l-g3)(l-g4)(l-a6) ' ^^^^ 



1 + 2g + 3g2 + 3g=^ + 2g^ 
(1 -g2)2(l -a3)2(i 

The complexity of the rational generators rapidly increases. For example 



^^"^ (l-g)(l-g2)(l-g3)(l-a4)(l-g6)(l-a9) ^ ^ 

has a numerator 

iVg (g) = 1 + 2g + 3g2 + 4g^ + 6g^ + 6g^ + 6g^ + la! + 4g^ + 5g^ + 4g^° + 2g" + 2g^2 - g^^ + . 

(16) 

As 72 ^ oo, a splendid behaviour of Dsk+2n,k follows directly from @: 

k\ lim ni-'^Z^afc+sn.fc = k lim(l - g)^Dfc(a) = (3 - V^)'^ + (3 + 73)"'= , (17) 

requiring in (|15|) that A^6(l) = 52, which mental arithmetic shows to agree with (|1^). 

Already a spectacular success emerges: the generator (|12[) is precisely that conjectured 
for depth-3 MZVs in [|], though no term of order occurs in the Ansatz (|^). Thus 
the simple-minded input of (|TU|,|rTp predicts (|T2p. The obvious question remains: how 
does the conjecture fare at greater depths? For A; > 3, data prior to BK was extremely 
scanty. Personal communication from Don Zagier, reporting computations performed by 
him and by Dror Bar-Natan, were consistent with the first 5 terms predicted by (0) at 
depth 4. Warmed by this success, we here extend it, by very large-scale computation, 
to the first 16 terms. The method involves shuffle identities p, |2^, outlined in 



the next section, which provide rigorous upper bounds on irreducibles at specific depths 
and weights. Confirming that 0^2,4 < HI involved 1.7 petashuffles. Having thus tuned 
and tested the code, we pushed on, to even larger systems of identities, at depth 5, where 
the first 7 terms of (0) have now been confirmed. The challenge of testing a significant 
number of terms in (pjf ) was judged imprudent to attempt, unaided. Fortunately help is 
at hand, from Jon Borwein, David Bradley and Roland Girgensohn, in a collaboration 
hosted by the Center for Experimental and Constructive Mathematics (CECM) at Simon 
Fraser University, which is a node of the Canadian High- Performance Computing Network. 



3. Depth-length and weight-length shuffles 



Depth-length shuffles were used in |25|, |26[. The idea is simple, in the extreme, 

though less trivial to notate. For simplicity, consider the product of a depth-1 Euler sum, 
whose argument string is merely {si}, and an Euler sum of depth k — 1, with argument 
string {s2,...,Sfc}. In the case of alternating Euler sums, signs may be included, as 
bars |1], Q or minus signs 0, in these strings. Then the product clearly entails the 
k Euler sums in which Si is inserted, in all possible ways, in the other string. The 
depth-length strings 

{Si, $2, ■ ■ ■ , Sk} , {S2, Si, . . . , Sk} , . . . , {S2, . . . , Si, Sk} , {^2, . . . , Sfc, Si} , 
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result, each with unit coefficient, together with sums of lesser depth, which do not concern 
the present analysis of reducibility. Generalization p8| to the product of a string of depth 
r with one of depth /c — r is as might be expected: take all shuffies of the two depth- 



length strings that preserve the order of each. Combined with the reducibility |^ of 

(^(argument string of length k) + (—1)'^^ (reversed string) 

such depth-length shuffles appear to exhaust that which can be concluded merely from the 
existence of a nested fc-fold sum; all else must take account of the form of the summand. 

In addition to such depth-length shuffle identities (called permutation identities in 0, 



24|, |25|) there are weight-length shuffle identities (corresponding to the partial- fraction 



identities of [0, ^ |25[). These follow from the existence of iterated- integral p2[ rep- 



resentations for Euler sums. In the case of an MZV of weight n, a representation has 
been given ^ |22|, 23| in terms of an n-fold iterated integral of the one-forms dx/x and 



dx/{l—x). The extension to alternating sums was given in [^]: one has merely to include 



the one-form dx/{l + x), which adds a third character to the weight-length alphabet. 

Weight-length shuffles are then expressions of the ring structure of alternating 
Euler sums, and their restriction to MZVs: 

C(string]^)^(string2) = ^ C (weight-length shuffled string) , (18) 

shuffles 

where each shuffled string preserves the order of each of the two constituent weight-length 
strings, in its iterated-integral representation. The product of strings of weights rii and 
n2 entails {^^^^^^ weight-length shuffles, each of which has a weight and a depth that is 
the sum of those in the product, and many of which may occur many times. 



4. Computational testing 



For double sums, (0) was the input to (H). However the correctness of that input is 
still, according to the author's understanding of p|, E4|, a conjecture. Thus it was tested. 



analytically, using REDUCE ||29[, to weight 44, which conffimed that the conjectured 
enumeration satisfied the identities of [p3], and required no other. 



For triple sums, (|T^ was a spectacular output of (^. However the correctness of that 
output is still, according to the author's understanding of 0, |2^, a conjecture. Thus 
it was tested, analytically, also using REDUCE, to weight 37, which conffimed that the 



conjectured enumeration satisfied the identities of |25|, and required no other 



The analytical results so obtained were consistent with the completeness and min- 
imality of the concrete bases proposed in Minimality is more or less unprovable, 
since nothing, in principle, forbids further reductions, beyond those entailed by shuffies. 
However, the lattice algorithm PSLQ |^ resolutely failed to find anything new, wherever 
we had the CPUtime and patience to probe. More such tests may now be made with 
PSLQ, since its author, David Bailey, has very recently implemented Richard Crandall's 



fast algorithm |31] for computing Euler sums, including alternations of sign, and achiev- 
ing a precision of 3200 digits at depth 5, and beyond. Conffimation has been obtained of 
instances of conjectures made in and in the course of the present work. Results may 
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shortly appear in the high-performance computing hterature. To date, none of the many 
numerical discoveries, made with PSLQ, erodes any conjecture in [jl], H, p6| . 

The strategy adopted to test (0) at depths 4 and 5, with high computational effi- 
ciency, was to pre-process all depth-length shuffles algebraically, using REDUCE, whose 
output was then converted to FORTRAN code. This was translated, by David Bailey's 
TRANSMP code, into multiple precision calls of MPFUN routines, which solved 



all the weight-length shuffles, modulo terms of lesser depth, enabling reducibility analysis 
at a numerical precision which provided overwhelming evidence that all known relations, 
and no others, were satisfied. The reader is (mercifully) spared programming details, 
which merely translate the transparent idea, of exhausting all shuffie identities, into com- 
putational practice. S/he is assured that the probability of misidentification of numerical 
solutions of the identities was less - often many orders of magnitude less - than 10~^°. 
The results that follow are hence not rigorously proven, though we discount the slen- 
der possibility that they fail to deliver the reductions that would have been achieved by 
computer algebra, with integer arithmetic, had one many gigabytes of core memory, and 
CPUyears of processing time, to expend on the huge integers that would be generated. 

A measure of what was involved is provided by the following statistics. At depth 
k = 4: and weight n = 42, the ring structure (|18]) generates 1720620718074180 shuffies 
(1.7 petashuffies) each of which contributes unity to an element of a 16000 x 10660 matrix 
of constraints, whose rank deficiency was found to be 111. At depth k = 5 and weight 
n = 27, the rank deficiency of a 29900 x 14950 matrix was found to be 36. 

In total, we confirmed the emboldened rank deficiencies in the integer sequence 

1, 1, 3, 5, 7, 11, 16, 20, 27, 35, 43, 54, 66, 78, 94, 111, 128, 150, 173, 196, 224, 254, 284, . . . 

(19) 



generated by ([T3|) , and those in the sequence 



1, 2, 5, 9, 15, 23, 36, 50, 71, 96, 127, 165, 213, 266, 333, 409, 498, 600, 720, 851, 1005, . . . 

(20) 

generated by (|l|). Neil Sloane has designated (|19|,|0D as A019449 and A019450 in the 
on-line encyclopedia of integer sequences. 

For reasons that may be appreciated, the CPUtime required for the above was consid- 
erable: about 20 CPUdays on a 256MB AlphaStation 600 5/333. For that reason, testing 
of (^ at depth 6 awaits input from colleagues at CECM p8|. Nor is this an idle pursuit. 



since at depth 6 there appear to be new |2^ features to solving the depth- and weight 



length shuffle identities. If (|]) is not correct, then depth 6 is the place to expect a failure. 
If a non-trivial number of terms in the amazing numerator (0) is, as hoped, confirmed, 
then the chances of failure at depths k > 7 appear remote, given the essentially depth-2 
nature of the input. 



5. Euler sums: a field theorist's sample bestiary 



Lest Euler sums and MZVs appear a diversion from serious-minded field theory, we briefiy 
recall some of their appearances in pQFT [|T], |^, |^, ^ |Tl], |1^, 0, |1^ |3^ , often in 
advance of their systematic investigation by the maths community. To this traffic, from 
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field theory to number theory, Dirk Kreimer has added an equally noble trade route: from 
field theory to knot theory 

Depth-1 Every pQFT practitioner knows that the odd zetas, C(2^ + 1), are every- 
where dense in multiloop QCD and QED. In counterterms, they correspond to 2-braid 
torus knots [16|. Here one beast is gloriously absent |l[ from pQFT: no sub divergence- 
free four-dimensional counterterm can spawn tt^, for which there is no knot. 

Depth-2 A decade ago IQ, before MZVs attracted apparent interest in the maths 



community, the author failed to reduce a depth-2 weight-8 Euler sum, in the e-expansion of 
dressed two-loop diagrams. The persistence ||15| of this state of affairs led to the prediction 
that it would surface in 6-loop counterterms. It does |0. It is the first depth-2 irreducible 
MZV, enumerated by 2 = 1. Even earlier [^, the St Petersburg group were unable 
to reduce e-expansions of critical exponents to depth-1. It is now known why: an infinite 
series [|lT| of depth-2 irreducibles occurs, starting at weight 8. The corresponding knots 
start at 8 crossings, with the 3-braid torus knot 819 = (4,3) p. Very recently p!7| , 
Anatoly Kotikov and the author have shown that alternating Euler sums result from the 
e-expansion of critical exponents, in = 3 — 25 dimensions. The first such beast is ({3, 1), 
which occurs with a multiple of C(4) that is process dependent. Thus odd- dimensional 
counterterms mimic even-dimensional massive diagrams, for it was observed in that 
two of the most important 3-loop results in pQFT - for the electron anomaly |3^ and the 
p-parameter |34| - have precisely that structure in their weight-4 terms. 



Depth-3 Likewise, there was an obstacle to the reduction of triple sums of weight 



11, obtained by cutting vacuum diagrams [T^. The irreducible occurs in the 7-loop beta- 
function of 0^-theory §]. As at depth 2 [Q, mathematicians later encountered 0, ^ 
this phenomenon. The associated knot is the sole positive 11-crossing 4-braid 0. Its 
uniqueness is required by -Dn^s = 1, and is confirmed by the knot enumeration of BK. 

Depth-4 Perhaps the most far-reaching gift from pQFT to number theory occurred 
at depth 4 and weight 12 0, [11|. Despite untiring effort, BK kept encountering an 



apparent mismatch between field, knot and (the then current state of) number theory. 



The (indubitably correct) result of |2^, that -Di2,2 < 1, jarred with the existence of a 
pair 1^ of 12-crossing 3-braid knots in the counterterms [|TT| of pQFT. The resolution 
was dramatic: the second knot corresponds to C(4,4, 2,2), which is pushed down |^ to 
C(9, 3), in the simpler enumeration (0), giving Mi2,2 = -Di2,2 + -Di2,4 = 1 + 1 = 2. Without 
this stimulus from pQFT, the author would not have discovered 0] such pushdowns, and 
the simplicity of (^ might have remained hidden. Instead of seeking a simple answer 
to the (now seen to be) difficult question Q3, BK backtracked to Q2, which has the 
simplest answer of the three. Then it was possible to return to Q3, input a minimal 
amount of information, interpreted in the light of Richard Feynman |jl|, [ll| and Vaughan 
Jones [0, observe the success ^ of (|^) in giving (0), and arrive, via (|T3|,|T^, at the 
predicted sequences of deficiencies in (|^,|^), tested here, non-trivially, at a cost of several 
petashuffles. 

Such are the fruits of computational field theory. 
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5. Conclusion and questions 



The conclusion is easy to state: the conjectured enumeration (||) emerged in the course 
of calculations in field theory [|I], ^, |1T|, is illuminated by knot theory 0, ||, 01; and is 
shown by p^, and by the emboldened results in ([I9|j20|) , to be in great shape, at 

all depths k < 6, however striking that may appear from the point of view of number 
theory g |6|, 1^. 



Further questions are easier to pose than to answer: 



Q4 Why? 

Q5 Is (|16|) correct, at depth 6? 

Q6 Is there an alternative to the CPU-intensive testing procedures, adopted here? 



These, and related issues, are receiving due attention B O, IT? . 
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